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Abstract 

We analyze the influence of pre-stress on the propagation of in- 
terfacial waves along the boundary of an incompressible hyperelastic 
half-space that is in contact with a viscous fluid extending to infinity 
in the adjoining half-space. 

One aim is to derive rigorously the incremental boundary condi- 
tions at the interface; this derivation is delicate because of the in- 
terplay between the Lagrangian and the Eulerian descriptions but is 
crucial for numerous problems concerned with the interaction between 
a compliant wall and a viscous fluid. A second aim of this work is to 
model the ultrasonic waves used in the assessment of aortic aneurysms, 
and here we find that for this purpose the half-space idealization is 
justified at high frequencies. A third goal is to shed some light on the 
stability behaviour in compression of the solid half-space, as compared 
with the situation in the absence of fluid; we find that the usual tech- 
nique of seeking standing waves solutions is not appropriate when the 
half-space is in contact with a fluid; in fact, a correct analysis reveals 
that the presence of a viscous fluid makes a compressed neo-Hookean 
half-space slightly more stable. 

For a wave travelling in a direction of principal strain, we obtain 
results for the case of a general (incompressible isotropic) strain-energy 
function. For a wave travelling parallel to the interface and in an 
arbitrary direction in a plane of principal strain, we specialize the 
analysis to the neo-Hookean strain-energy function. 
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1 Introduction 



Seismic records show that underground rocks and ocean beds are subject to 
stress and strain and that surrounding fluids are viscous and under high pres- 
sures. Chnical ultrasonic measurements indicate that arteries can undergo 
large strains in service and are most sensitive to changes in blood pressure. 
Many moving and vibrating parts of automotive devices are made of loaded 
elastomers in contact with highly viscous fluids. These are a few examples of 
situations where it is crucial to model and understand the motions and the 
stability of the interface between a deformed elastic solid and a viscous fluid. 
Yet only a handful of studies can be found on the subject, especially when 
compared with the abundant literature on waves at the interface between an 
elastic solid and an inviscid fluid, which goes from the pioneering works of 
Galbrun yj, Cagniard [2J, Scholte [3], and Biot [1] to the definitive treatment 
of the acoustoelastic effect by Sinha et al. [5]; see also Poiree [6] and Degt- 
yar and Rohklin [7]. Waves at the interface between a viscous fluid and an 
undeformed isotropic elastic solid were examined by Vol'kenshtein and Levin 
[8j and the corresponding problem for an anisotropic elastic solid by Wu and 
Wu [9] . To the best of our knowledge, only Bagno, Guz, and their co-workers 
have studied the title problem (see, for example, [10], |II])- Their analytical 
treatment is, however, quite succinct and we therefore aim to shed new light 
on the problem by re-examining it on the basis of recent developments in the 
theory of small-amplitude waves, linearized in the neighbourhood of a finite, 
static, homogeneous deformation. 

It turns out that one of the trickiest aspects of the study is the derivation 
of proper incremental boundary conditions at the interface because these are 
usually written in terms of the nominal stress in a deformed solid (Lagrangian 
formulation. Section 2) and in terms of the Cauchy stress in a fluid (Eulerian 
formulation. Section 3). These equations are combined in an appropriate 
way for a general interface in Section 4. We then specialize the analysis to 
principal wave propagation for an arbitrary (incompressible, isotropic) strain- 
energy function in Section 5. In the course of the analysis in Section 5, by way 
of application of the theory, we show that in respect of an abdominal aortic 
aneurysm it is appropriate to neglect the curvature and finite thickness for 
ultrasonic waves (10 MHz), i.e. to treat the aneurysm locally as a half-space. 
It also shown that it is not appropriate to use waves with a real frequency 
to study the stabihty of compressed solids in contact with a viscous fluid. 
Finally, in Section 6, by specializing to the neo-Hookean solid, we consider 
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the propagation of non-principal waves for both tension and compression of 
the half-space in order to illustrate the influence of the fluid. 



2 Basic equations for the solid 

For the solid material we denote by F the deformation gradient relating 
the stress-free reference configuration, denoted Eq, to the finitely deforming 
configuration, denoted E. This has the form F = Grada:;, where x = x(X, t) 
is the position vector in E at time t of a material point located at X in ;Bo, 
X is the deformation mapping, and Grad is the gradient operator relative to 
Eo- 

We consider the material to be elastic with a strain-energy function, de- 
fined per unit volume, denoted hj W = W{F). Furthermore, we restrict 
attention to incompressible materials so that the constraint 

detF = l (2.1) 

is in force. The nominal stress tensor, here denoted by S, and the Cauchy 
stress tensor cr are then given by 

dW 1 dW 

S = ^lL- pF-\ a = F^ - pi, (2.2) 
OF ^ ' OF ^ ' ^ ^ 

where p is a Lagrange multiplier associated with the constraint (12.11) and I 
is the identity tensor. 

The equation of motion is 

Div5 = px, (2.3) 

where Div is the divergence operator relative to Eq, p is the mass density of 
the material, and a superposed dot signifies the material time derivative. 

Next, we consider a small motion superimposed on the finite deformation. 
Let u{x,t) be the displacement vector relative to E and v{x,t) = ii the 
associated particle velocity (the material time derivative of u). Then, on 
taking the increment of equation (12. 3p and thereafter changing the reference 
configuration from Eq to E, we obtain 

divs = pu = pv, (2.4) 
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where s is the increment in S (referred to B) and div the divergence operator 
relative to B. 

The (hnearized) incremental version of the constitutive relation (12.21) is 
written 

s = A(gradtt) + p(grad'u) — p7, (2.5) 

where A is a fourth-order tensor of elastic moduli, grad is the gradient 
operator relative to B, and p is the increment in p. In component form, this 
is written 

Sij = Aijkiui^k + pui,j - p6ij, (2.6) 

where denotes d/dxk and 6ij is the Kronecker delta. In terms of W the 
components of A. are given by 

A,ki = F,pFk,^^r^^- (2.7) 

For details of these derivations (in a slightly different notation) we refer to 
Dowaikh and Ogden [12] . 

We now consider the material to be isotropic, so that W = W{Xi, A2, A3) 
is a symmetric function of the principal stretches, Ai,A2,A3 (the positive 
square roots of the principal values of FF"^ , where signifies the transpose), 
subject to the constraint 

A1A2A3 = 1, (2.8) 

which follows from (12. ip . Then, on noting that for an isotropic material cr is 
coaxial with FF^ and specializing equation (I2.2p 9. we obtain the principal 
Cauchy stresses (see, for example, Ogden [13] ) 

(Tj = —p + \iWi, = 1, 2, 3 (no sum over i), (2.9) 

where Wi = dW/d\i, i = 1, 2, 3. 

When referred to the same principal axes, the only non-zero components 
of A. are 

Aiij, = XiXjWij, t,j e {1, 2, 3}, (2.10) 
Am = ^^^^^f^^^l ^,jG{l,2,3},.^j, (2.11) 

A-ijji Ajiij A-ijij AjPVi, Z, J G '[1, 2, 3]-, i ^ J, (2.12) 
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where Wij = d'^W/dXidXj. For subsequent convenience, we adopt the nota- 
tions defined by 

7., = (W,-A,4y,)A,V(A,2-A,2), 

= (X-Wu + A|iy,,)/2 - X^XjW^, + (X^W, - X,W^)X,X,/iX^ - A|), (2.13) 

noting that 7jiAi = lij^] and Pji = Pij. 

2.1 The pre-stressed elastic half-space 

We now consider B to be independent of time and to correspond to a pure 
homogeneous strain of a half-space defined by X2 ^ 0. The half-space is 
maintained in this configuration so that its boundary X2 = is a principal 
plane of strain. We denote by Xi and the other two principal directions of 
strain and by Ai, A2, A3 the principal stretches in the Xi, X2, X3 directions, 
respectively. The corresponding principal Cauchy stresses are then as given 
by (12.91) . In particular, the boundary 0:2 = is subject to a normal stress (T2 
and, after elimination of p, the other two principal Cauchy stresses are then 
given by 

ai = a2 + XiWi - X2W2, as = a2 + X3W3 - X2W2. (2.14) 

We are interested in the propagation of incremental (small amplitude) 
acoustic waves along the boundary plane X2 = 0, in a direction making an 
angle 6 with the principal direction xi. The incremental velocity and nominal 
stress fields v and s are then considered as superimposed on this finite static 
configuration. We examine inhomogeneous time-harmonic plane waves of the 
form 

{v, s}(xi, X2, X3, t) = {v{x2), -(A;/cu)s(x2)}e''=(^«^^+^«^^)e-''^*, (2.15) 

where we have introduced the notations ce = cos 9, se = sin 9, k is the wave 
number, uj is the wave frequency, and v, s are functions of X2 only, such that 

v{oo) = 0, s(oo) = 0. (2.16) 

Using the results of Destrade et al. [14] (see also Chadwick pE]), we find that 
the incremental equations of motion can be written as a first-order differential 
system of six equations, namely 

^'{x2)=ikN^{x2), (2.17) 
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where the notation ^ is defined by 



€ = [Vi,V2,V3,S21,S22,S23\ 

and tlie 6x6 matrix N has the block structure 



N 



A^i N2 
N3 + PI 



(2.18) 



(2.19) 



in which the 3x3 matrices Ni, N2, N3 are real and their components depend 
on the material parameters 'jij and /3ij given in fl2.13p . and the notation 
p = pu^/k"^ has been introduced. Here / represents the 3x3 identity matrix. 
Explicitly, -Ni, N2, -N3 are 



(2.20) 



ce(72i 


- 0-2) /721 


0" 




1/721 







1] —K 







se 












z/ 


se(723 


- 0"2)/723 


0_ 







1/723. 




— K, p 



respectively, where 
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Se73i, 



T] = 2cq{(3i2 + 721 - 0-2) 
^ = Ce[7i2 - (721 - 0-2)7721] + sl\j32 - (723 - 0-2)7723], 
P = C0713 + 2sg(/923 + 723 - 0-2), 
K = ceSeiPis - P12 - P23 - 721 - 723 + 2cr2). 



(2.21) 



Equation fl2.17p provides the general expression for the equations of motion, 
for arbitrary 6 and W. 

Now, in seeking a decaying partial-mode solution of the form 



-ksxi 



C, 3ft(M > 0, 



(2.22) 



where ^ is a constant vector and s an unknown scalar, we arrive at the eigen- 
value problem N(^ = is^. In general, the associated propagation condition, 
det{N—isI) = 0, is a cubic in [12], where now / is the 6x6 identity matrix. 
Its analytical resolution is too cumbersome to be of practical interest, and so 
we specialize the general equations to the following, simpler, situations: (i) 
principal wave propagation [6 = 0) for arbitrary W; (ii) non-principal wave 
propagation {6 7^ 0) for the neo-Hookean material, for which 



W = C{Xl + Xj + Xl - 3)/2 



(2.23) 
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where C > is a constant (the shear modulus of the material in the reference 
configuration) . 

In Case (i) , the equations of motion decouple the system 



'23. 



ik 



1/721 

1/723 



^3 
^23 



(2.24) 



(for which the trivial solution may be chosen) from a system of four differen- 
tial equations for vi, V2, S21, S22- Hence, in this case, the wave is elliptically 
polarized, in the {xi,X2) plane. The corresponding propagation condition is 
a quadratic in s^, which can be solved explicitly. 

In Case (ii), we also find that the wave is two-partial, polarized in the 
plane containing the directions of propagation and attenuation (the saggital 
plane); there, the corresponding propagation condition involves the product 
of the factor (s^ — 1) and a term linear in s^, which simplifies the analysis. 

Before embarking on the details of these cases, we complete the descrip- 
tion of the boundary- value problem by considering the behaviour of the wave 
in the fluid in the half-space X2 <0. 



3 The fluid half-space 

Adjoining the deformed solid half-space is a half-space X2 ^ filled with an 
incompressible viscous Newtonian fiuid, for which all mechanical fields are 
denoted by a superscripted asterisk. In the static state the fiuid is subject 
only to a hydrostatic stress cr* — —P*I, and by continuity of traction across 
the boundary X2 = we must have 

- P* = (T2. (3.1) 

The constitutive law for the fiuid associated with the motion is then written 
in terms of a superimposed Cauchy stress tensor, denoted here by s* and 
given by 

s* = -p*I + 2^*D\ tri?* = 0, (3.2) 
where 11* is the viscosity of the fiuid, 

Z)* = i[gradt;* + (grad7;*)^], (3.3) 

V* is the fiuid velocity, and p* — p*{x, t). 
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We seek inhomogeneous waves with the same structure as in the sohd, 
that is 

{v*, s*}{x,, X2, X3, t) = {v*{x2), -(A;/a;)s*(x2)}e^'=(^«^i+^«^3)e-''^*, (3.4) 

where v*, s* are functions of X2 only, such that 

v*{-oo) = 0, s*(-oo) = 0. (3.5) 

We find that the equations of motion, divs* = p*v* (where p* is the mass 
density of the fluid), linearized in v*, can be cast as 



e'{x2) = ikN*e{x2), 



where 



1^* r^* ^* ^* ^* ^* ^* iT" 

S — '^2) '^3' *21) *22) *23J ' 

and the constant complex matrix N* has the block structure 

Ni n; 



N* 



Nl + p*I AT* 



(3.6) 
(3.7) 

(3.8) 



"0 


Ce 


0' 




"1 0" 




'44 + si 





2>cese 


ce 





Se 





















se 







1 




Zcese 








N2, being real symmetric matrices, and the notation p* = p*uj'^ /k"^ 
has been adopted. If we write fi* = p*u!, then, respectively, —Nl, —ip,*Nl, 
-iNl/fi* arc 



(3.9) 



Again, when we seek a decaying partial-mode solution, this time in the 
form 

^*(X2) = e^''^X*, ^(ks*) > 0, (3.10) 

where is a constant vector and s* an unknown scalar, we end up with 
an eigenvalue problem, here N*(^* — — is*^*. The associated propagation 
condition is det(iV* + is* J) = 0, which here simplifies to 



with roots 



(s*'-l)(s*'-l + ie)^ = 0. 



± 1, ±-\/l — ie (repeated). 



(3.11) 
(3.12) 
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where e = p* / ft* = p*u!/{fj,*k'^). The roots are independent of 9, as expected, 
because the fluid is isotropic. Corresponding to each of the four roots, there 
are four eigenvectors and therefore potentially four partial-modes. However, 
two of these must be discarded since their amplitudes do not decay with 
distance from the interface X2 — 0. The two remaining modes form the basis 
for the general solution of the equations of motion that is needed for matching 
with the two-partial wave in the solid. 

We now give the general boundary conditions at the deformed solid/ viscous 
fluid interface. 

4 The interface 

In order to match the incremental tractions across the boundary it is nec- 
essary to work in terms of the Cauchy stress since the nominal stress is not 
deflned inside the fluid. Towards this end we flrst calculate the incremental 
traction in the solid in terms of the Cauchy stress. Continuity of traction 
requires 

S^NdA = anda = a*nda, (4.1) 

where dA and da are the area elements in Bq and B, respectively. Taking 
the increment of this and updating the reference configuration to B yields 

s^nda = anda + anda = s*nda + (7*nda, (4.2) 

where a superposed tilde indicates an increment. Note that, after updating, 
F = I and S = cr in the configuration B. 

Now, according to Nanson's formula (applied to the boundary of the 
solid), we have nda — F~^NdA, from which it follows, again after updating, 

that 

ndd — — (gradw)^nda. (4.3) 

Hence, the incremental traction continuity condition can be written 

s^n = [a- — cr(gradw)^]n = [cr* — (T*(gradw)'^]n, (4.4) 

and we recall that cr* = —P*I. 

Since n is in the X2 direction for the considered half-space we may write 
the continuity condition in component form as 

S2i = + P*U2,i, « = 1, 2, 3, on X2 = 0. (4.5) 
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Additionally, the velocity must be continuous, i.e. 

V* = Vi, i = 1, 2, 3, on X2 = 0. 



(4.6) 



In terms of the functions v{x2), s{x2) and their counterparts in the fluid, 
the boundary conditions become 

v*{0)=v,{0), t = 1,2,3, (4.7) 

and, noting that Vi = —iuui, 

sl^{Q) + ceP*V2{0) = 52i(0), sUO) + seP*V2iO) = §23(0), (4.8) 

and 

- ^P%{0) = §22(0). (4.9) 

5 Principal waves: no restriction on W 
5.1 General solution in the solid 

Here we take ^ = and place no restriction on the form of W. When 6 = 0, 
the in-plane mechanical fields in the solid satisfy the equations of motion 
^' = ikN^, where now |(x2) = [vi,V2, hi, -§22]'^ and 



N 






-1+^2 


1/721 





-1 











— 1] 








-1 





p- z/ 


-1 + 0^2 






(5.1) 



with (T2 = 0-2/121 and t] and u now reduced to 

T] = 2[Pi2 + 721(1 - CT2)], «^ = 712 + 721(1 - 0^2)^ 
Here the propagation condition is a quadratic in [12] and given by 



where 



2(3 



2(3i2 - P 
721 



a 



0, 



7i2 - P 
721 



(5.2) 
(5.3) 
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We recall that p = pu'^/k'^. Formally, the roots of the quartic fl5.2p are 

± s/JV^fW^, ±^Jl3-^W^- (5.4) 

We pause the analysis to highlight a particular feature of the present 
interface waves. Because the fluid is viscous, the wave number k is complex 
and so, therefore, are j3 and a. It follows that, in contrast to the purely 
elastic case [T7|, it is not clear a priori which two of these four roots are such 
that the decay condition Eq. fl2.22p 9 is satisfied. Let Si and S2 be two such 
roots. We note first that 

sl + sl = 2(3, sjsl = a\ (5.5) 

and then, depending on which value in Eq. (15.41) they correspond to, one of 
the following four possibilities may arise: 



S1S2 = a. 


Sl + S2 = 


S1S2 = a, 


Sl + S2 = 


S1S2 = -a, 


Sl + S2 = 


S1S2 = -a, 


Sl + S2 = 



V2 {(3 + a), 
-^2{(3 + a), 
V2(/5-«), 

-^/2{(3-a). (5.6) 

Now, the eigenvectors and associated with si and S2 are columns 
of the matrices adjoint to N — isi J and N — 182!, respectively. Taking, for 
instance, the fourth column gives 

C = [is., -1, -721(1 + - aa), -i72is.(l + 2/? - - ^2)]^, 2 = 1,2. (5.7) 

Finally, the general solution in the solid is 

^(X2) = Aie-'^i^^C' + Ase-'^^^^^C', (5.8) 
where Ai and A2 are constants. 

5.2 General solution in the fluid 

In the fluid, a two-partial solution is required for adequate matching at the 
interface. Taking 6* = in Section 3, we find that the in-plane equations of 
motion are £*' = ikN*^*, where now ^*(x2) = [vl,V2, §21, §22]'^ and 

-1 i/fx* ' 
-1 
Aip* + p* 0-1 

p* -1 



N* 



(5.9) 
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We assume that the wave propagates and is attenuated in the direction 
Xi ^ 0. From Eq. fl3.4p with 6 = 0, we see that these assumptions lead to 



^{k) > 0, Q{k) > 0. 



(5.10) 



On the other hand, the wave is also attenuated with distance from the inter- 
face; it follows that we can discard the root —1 from Eq. (13.121) and retain 
the root +1. The choice to be made for the remaining two roots in Eq. (13.121) 
is not so clear cut and for the time being we call s* the suitable root; hence 
s* is such that 



s" = 1 — ie, ^{ks 
Finally, the general solution in the fluid is 



(5.11) 



where and A2 are constants 



1 
1 

-2i/i* 

-/i*(l + S*2) 



+ A;e^'*^^ 



1 

-i/i*(l + S*2) 

-2fi*s* 



(5.12) 



5.3 Dispersion equation for the interface wave 

When we specialize the general boundary conditions (14.71) . (14. 8 p and (14. 9 p to 
the present context, we find a linear homogeneous system of four equations 
for the four unknowns Ai, A2, Al, A2, for which the associated determinant 
must be zero. After some manipulations, using Eq. (15.50 . we find that 



-Si -S2 1 S* 

1111 
721(1 + S?) 721(1 + si) -2i/i* -i;i*(l + s^2^ 

72lSl(l + s2) 72152(1 + 5?) ifi*{l + S*^) 2ifl*S* 



0. (5.13) 



We see at once that the normal load a2 does not appear explicitly in this 
equation. This feature highlights a major difference between a wave at the 
interface of a loaded solid half-space and a vacuum (Rayleigh wave) and 
a wave at the interface of a loaded solid half-space and a viscous fluid, as 
considered here; Chadwick and Jarvis [18] also noted this peculiarity for 
waves at the interface of two loaded solid half-spaces (Stoneley wave). Of 
course, a2 still plays an important role, in particular in the determination of 
the pre-stretch ratios and of the amplitudes of the tractions in the solid. 
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Once the determinant is expanded and the factors (si — S2)(l — s*) are 
removed, we end up with the exact dispersion equation, 



7i2 - p , 2(3i2 + 2721 - p 

\ S1S2 - 1 

721 721 

- i72l/i*[2(l - SiS2)(l - S*) + (Si + S2){1 + S*){S* + S1S2)] 

+ [ie + (ie - 4) s*] = 0. (5.14) 

This equation is fully explicit because the terms si + S2, S1S2, and s* are 
given by Eqs. ( I5.6P and Eq. (15. lip . Of course, as noted earlier, there are 
several possibilities for these terms, which generate in total eight different 
dispersion equations. In each case the resulting root(s) in k must be checked 
for validity against the propagation and decay conditions (I2.22p 9. f l5.10p . and 
f l5.1ip 9. which we summarize here as 

^{k) > 0, '^{k) > 0, ^{ks*) > 0, ^{ksi) > 0, ^{ks2) > 0. (5.15) 

Notice that in the special case of solids whose strain-energy function W is 
such that 2Pi2 = 712 + 721, which includes the neo-Hookean solid fl2.23p . the 
biquadratic (15.20 factorizes as (s^ — l)(s^ — a^) = 0. Hence si = 1, S2 = ±a, 
and the exact secular equation (15.140 simplifies accordingly, leading this time 
to four different exact dispersion equations. 



5.4 Application: modelling of intravascular ultrasound 

In recent years, intravascular ultrasound (IVUS) has proved to be a most 
promising tool of investigation for measuring and assessing abdominal aortic 
aneurysms. Its accuracy is as good as that of computed tomography (CT) 
scans and it has obvious non-radiative advantages [12], [20], [21] ■ We now 
apply the results of this section to an IVUS context. 

First we recall that medical ultrasound imaging devices operate in the 
1-10 MHz range; accordingly we take u = 10^ Hz. We argue, and will check 
a posteriori, that at such high frequency the wavelength is small compared 
with the radius and thickness of an artery so that, as far as the propagation 
of localized waves is concerned, an aortic aneurysm can be modelled as a 
half-space. Here we take the axis of the artery along the xi direction and 
consider that the half-space 0:2 ^ is filled with blood. 

For the solid, we use a strain-energy function devised by Raghavan and 
Vorp [22J to fit experimental data collected on uniaxial tension tests of aortic 
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aneurysms, namely 

W = Ci{\\ + + - 3) + C2{\1 + + 3)2, (5.16) 

where, typically, Ci = 0.175 MPa, C2 = 1.9 MPa. We assume that the 
aneurysm corresponds to a region where the tissue undergoes a strain of 
20%, and that the end-systolic blood pressure is high, 150mmHg (= 20kPa) 
say, so we take (T2 = 20 kPa. We consider that the tissue is free to expand or 
contract in the X3 direction with = 0. From (I2.14p we find the remaining 
stresses and strains. Summarizing, we have 

Ai = 1.2, A3 - 0.908, ai~ 71.45 kPa, a2 = 20kPa, (T3 = 0, (5.17) 

with A2 calculated from the incompressibility condition. For the mass density, 
wc take p = 1000 kg/m^. For the blood, we use typical values of viscosity 
[23J and mass density |21j: fx* = 3.5 x lO'^Ns/m^, p* = 1050 kg/m^. 

Here, the frequency oo is fixed as a real quantity a priori. We then replace 
k everywhere hj k = ujS, where S = + iS~ is the (complex) slowness, 
and the only unknown in the dispersion equation fl5.14p . The propagation 
and decay conditions (15.151) are satisfied and may be written as 

^{S) > 0, 53(5) > 0, ^{Ss*) > 0, ^{Ssi) > 0, ^{Ss2) > 0. (5.18) 

We find that the only qualifying root is 

5 = 2.721 X 10^2 + 4.424 X 10"^i, (5.19) 

from which we deduce the phase speed v = 1/5'+, the damping factor 7 = 5~, 
and the wavelength Aq = 27r/(ti;5"'"): 

u = 36.75 m/s, 7 = 4.424 x lO'^m'S Aq = 23.1 /im. (5.20) 

Then we plot the depth profiles of the wave. Its amplitude is the real 
part of V in the solid and of v* in the fluid. Explicitly, 

3? {v{xi,X2,t)) = e"'^'^'' [3? {v{x2)) cosuj{xi/v - t)] . 

—Q{v{x2))smuj{xi/v — t)], (5.21) 

in the solid, and similarly for the fluid, for which v and v are replaced by their 
asterisk counterparts. Clearly, the particle velocity is elliptically polarized. 
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and the lengths of the elhpse semi-axes decay with distance away from the 
interface and also with increasing Xi. Figure 1 shows the variations of the 
normal (continuous curve) and tangential (dotted curve) velocity components 
in the fluid {x2 ^ 0) and in the solid (x2 ^ 0), normalized with respect 
to ^2(0), as functions of X2/A0. Note that the components are in phase 
quadrature. Note also, as is clear from the zooms shown in Figure 1, that 
the components are continuous across the interface, as expected, but their 
first derivatives are discontinuous. The wave is elliptically polarized near the 
interface; the major axis is normal to the interface and more than 12 times the 
length of the minor axis; the ellipse is described in a retrograde manner. This 
shape is carried through the depth of the solid whereas in the fluid, the wave 
becomes rapidly circularly polarized at a depth of about 0.06 wavelengths, 
and remains nearly so through the rest of the half-space. The localization is 
greater in the fluid than in the solid: the amplitude has almost vanished after 
one wavelength into the former and after five wavelengths in the latter. An 
aneurysm is typically 1 mm thick, which, with the numerical values used here, 
is more than 50 wavelengths; thus, the assumption of a semi-infinite solid is 
justified, as is the assumption of a flat interface (aneurysms are typically of 
diameter 5 cm). 

Finally, we note that if the fluid is absent then the corresponding surface 
wave would travel with speed 40.734 m/s; hence the viscous fluid not only 
dampens the wave but also slows it down noticeably. Of course, blood flows 
in an artery, and creates shear deformation and stress in the solid. The 
influence of wall shear stress on the waves will be treated elsewhere. 

5.5 Example: compressive stresses 

Here we consider the behaviour of the elastic half-space (in contact with the 
viscous fluid) when it is compressed in the xi direction (Ai < 1). In their 
pioneering works. Green and Zerna [25] and Biot [26], [2^ showed that when a 
highly elastic half-space with a free surface (no fluid contact) is compressed, 
a surface instability may develop. Focusing on neo-Hookean solids 02.231) 
they showed that the critical stretches are Acr = 0.666 for tangential equi- 
biaxial compression (Ai = A, A2 = A~^, A3 = A), Acr = 0.544 for plane strain 
compression (Ai = A, A2 = A^^, A3 = 1), and A^ = 0.444 for normal equi- 
biaxial compression (Ai = A, A2 = A~^/^, A3 = A^^/^). Theirs was a static 
stability analysis, later also included in a wider dynamical context by Flavin 
[28], Willson [2n], Chadwick and Jarvis [30], Dowaikh and Ogden [12], and 
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Figure 1: Depth profiles as functions of X2/X0 of tlie velocity components of 
the acoustic wave: X2 (normal) component - continuous curve; Xi (tangential) 
component - dotted curve. The zooms show continuity of these components 
and discontinuity of their derivatives. 



others. We now consider this problem. 

A localized small-amplitude wave propagates over the free surface of a 
deformed Mooney-Rivlin or neo-Hookean half-space with normalized squared 
speed pv'^/'ju = 1 — cto'^I/'^ii where Uo is the real root of cr^ + cr^ + 3(t — 1 = 
((Jo = 0.2956). Clearly, in the examples of plane strain and equi-biaxial strain 
above, the squared wave speed increases when A increases and decreases when 
A decreases. At A = 1 there is no pre-strain and pv'^/C = 0.9126, the value 
found by Lord Rayleigh [31] for isotropic incompressible linearly elastic solids. 
As A — > 00, the squared wave speed tends to the squared wave speed 712/p 
of a transverse bulk wave. As A decreases, there is a critical stretch Acr at 
which the squared speed is zero and below which f ^ < 0. Since the wave time 
dependence is of the form Q^^i^i-'"^^) ^ with A; > 0, it follows that the amplitude 
grows without bound in time when A < Acr, and that the surface becomes 
unstable (at least, in the linearized theory). It is therefore appropriate to ask 
what happens when the compressed half-space is in contact with a viscous 
fluid. 

Bagno and co-workers addressed this question in a series of articles (see 
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the review by Bagno and Guz [TT] and references therein) . For a neo-Hookean 
sohd they found that the wave speed drops to zero when A = 0.544 in plane 
strain compression and when A = 0.444 in normal equi-biaxial compression, 
that is at the same critical stretches as for surface (solid/ vacuum) instability, 
irrespective of the viscous fluid characteristics. On the other hand, Bagno 
and Guz [TU] find that the wave speed falls to zero at a critical stretch that 
does depend upon the viscosity of the fluid. To address this disparity, we 
now compute the speed of the interfacial wave when the half-space is in 
compression. 

In order to minimize the number of parameters, we use the neo-Hookean 
solid, with W given by (12.231) . We take a normal equi-biaxial pre-strain 
(Ai = A, A2 = A~^/^, A3 = A^^/^) and, following Bagno and Guz, we take the 
frequency uj to be real. A dimensional analysis of the resulting dispersion 
equation shows that (I5.14p now depends on just three non-dimensional pa- 
rameters: a measure of the pre-strain. A; a measure of the dynamic viscosity 
of the fluid compared with the shear modulus of the solid, ii*uj/C] and the 
ratio of the densities, p* / p. Once these quantities are specified, the dispersion 
equation may be solved for the non-dimensional complex unknown x defined 
by 

where S is the (complex) scalar slowness. The dispersion equation can now 
be solved numerically for x, and the interfacial wave speed, normalized with 
respect to the transverse bulk shear wave speed in the deformed solid, is 
c = 3fJ(l/x). 

For Figure 2(a), we fix p* / p at 1.0 and we take in turn p*uj/C = 0.2, 0.04, 
0.02. The first choice {p* / p = 1.0, p*uj/C = 0.2) is roughly that obtained for 
the blood/artery interface of the previous section with Ci = C and C2 = 0. 
We see clearly that as p*uj/C decreases the wave speed decreases towards 
zero as A tends to 0.444, the critical compression stretch for the solid/vacuum 
interface (the thick curve gives the solid/vacuum interface wave speed). We 
find that in the extension to moderate compression range, the solid/fluid 
interface wave speed is significantly lower than the solid/vacuum wave speed. 
In the strongly compressive range (as A 0.444), the speed plot dips towards 
zero, and dips further as p*uj/C decreases, without ever reaching that value] 
the plot gets close to the plot for the solid/vacuum interface wave speed 
as A decreases, but then crosses it, and the fluid/solid interface wave speed 
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Figure 2: Normalized wave speed for a neo-Hookean solid under normal equi- 
biaxial pre-strain and a viscous fluid: (a) p* / p = 1.0 and p*u}/C = 0.2, 0.04, 
0.02; (b) p*/p = 1.0, 0.4, 0.2 and p*u/C = 0.1. The thick curves are for an 
unloaded solid half-space. 

increases again. Note that we checked the validity of the solution at all 
compressive stretches using (15.181) . 

For Figure 2(b), we fixed p*uj/C at 1.0 and took p*/p = 0.1,0.05,0.01. 
In the extension to moderate compression range, the plot for the solid/fluid 
interface wave speed gets closer to that for the solid/vacuum interface wave 
speed as p*/p decreases. In the strongly compressive range, similar comments 
to those made for Figure 2(a) apply. 

The conclusion is that as both p*uj/C and p* / p tend to zero, the speed 
tends to zero when the stretch tends to the critical compression stretch of 
the solid/vacuum surface instability. This is to be expected because this 
double limit corresponds to the vanishing of the fluid. However, as em- 
phasized earlier, the speed never reaches the zero limit and the bifurcation 
criterion is therefore never met. For instance, the typical values used by 
Bagno and co-workers [H] put p*uj/C at about 0.0002 and p* / p as low as 
0.1, giving c smaller than 10~^, but not zero. In effect a localized damped 
wave exists for the whole compressive range, with speed values starting below 
the solid/vacuum interface wave speed in the moderate compressive range, 
reaching a minimum (above the solid/vacuum interface wave speed) as A ap- 
proaches 0.444, and then rising rapidly to infinity as A decreases below 0.444. 
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This result suggests that when a neo-Hookean half-space is in contact with 
a viscous fluid it becomes completely stable. This, however, is an incorrect 
deduction, because it is based on special motions, for which the frequency lo 
is assumed real. Other motions might be unstable, as is illustrated below. 

We now take the wave number k to be real {k > 0) and let the speed v 
be complex, 

> 0, v = v+ + iv-, (5.23) 

so that the motion is now proportional to e'^'' t^ik{xi-v+t) ^ Clearly in this 
case, the conditions for a stable, localized wave, propagating in the xi > 
direction, are 

^{v) ^ 0, '^{v) ^ 0, ^{s*) > 0, 3fJ(si) > 0, 3f?(s2) > 0. (5.24) 

(note that 3fJ(si) > is automatically satisfied in a neo-Hookean solid, be- 
cause Si = 1). 

Then, an analysis of the dispersion equation (15.141) . in the case of a 
neo-Hookean solid with normal equi-biaxial pre-strain, reveals three non- 
dimensional quantities: A, ii*k/ pC , and p* / p. Once they are specified, we 
solve the dispersion equation for the non-dimensional quantity x defined as 




(5.25) 



The interfacial wave speed, normalized with respect to the transverse bulk 
shear wave speed in the deformed solid, is c = 3fJ(x). 

For Figure 3, we take p* / p = 1.0, and p*k/^/pC = 0.2, 0.002 in turn, and 
we plot V against A in the compressive range. We find that at a compressive 
stretch close to the critical compressive stretch of surface stability for the 
solid/vacuum interface, the normalized speed c drops to zero. From the 
figure it is not clear that the values of A at this point are different for 0.2 
and 0.002, but the zoom in Figure 3 shows, however, that the value of the 
compressive stretch at which c = depends on the material parameters. 
For comparison, the curve corresponding to p*k/ y/pC = 1 is also shown. We 
emphasize that the situation at c = does not correspond to a static solution 
of the equations of motion (which would be impossible in the fluid), but to 
a non-propagating damped motion, proportional to e'^'^ t^ikxi^ Moreover, 
this situation does not correspond to an instability because at that point, 
and a bit beyond, the requirements f l5.24p still hold. For instance, in the case 
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Figure 3: Normalized wave speed for a neo-Hookean solid under normal 
equi-biaxial pre-strain and a viscous fluid with p* / p = 1.0 and ix*k/\/pC = 
0.2, 0.002, 1 (indicated by the arrows) for the compressive range 0.4 < A < 1. 
The zoom shows the details for p*k/y/pC = 0.2, 0.002 for 0.442 < A < 0.445 
range. The thick curve corresponds to the solid half-space with no fluid 
loading. 

p*/p = 1.0, p*k/y/pC = 0.2, the speed drops to zero at A ~ 0.44539 but there 
still exist non-propagating, localized motions for A > 0.42212. Beyond that 
stretch value, however, all solutions of the considered type grow unboundedly 
with time and/or with space, indicating instability. From this example we 
see that the viscous load stabilizes slightly the solid half-space, because it 
can now be compressed by an extra 2%, from 0.444 to 0.422. 

6 Non-principal wave: neo-Hookean solid 

In this section we consider that the solid half-space is a deformed neo- 
Hookean material, with strain-energy function given by Eq. ( 12.23^ . For this 
material. Flavin p8j noticed that Rayleigh waves are plane-polarized for any 
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direction of propagation in a principal plane; see [22], [23] • Chadwick and 
Jarvis [TBj showed the same result for Stoneley waves. Indeed, we now show 
that the saggital motion is always decoupled from the 'anti-saggital' motion, 
for any 6 and any type of (bulk or interfacial) inhomogeneous wave. 

First, we note that for neo-Hookean bodies, the parameters •jij and Pij 
have the simple expressions 



7., = CXl A, = C(A^ + A?)/2. 



(6.1) 



It follows that the matrices Ni, N2, N3, given by Eq. fl2.20p . are greatly 
simplified. Next, we recall that the direction of propagation and the normal 
to the interface define the saggital plane, with unit normal [—sq, 0, cq]^. Now 
consider the new unknown functions Wi, (^ = 1, 2, 3), defined by 



Wi 



where 



Cg Se 

1 

-Se Cg 



(6.2) 



(6.3) 



Some algebraic manipulations reveal that the equations of motion (12.1 7p . 
written for Wi, t2i, decouple the 'anti-saggital' motion [1^3,^23] from its sag- 
gital counterpart. For the latter we find 



[w[, w'2, t'ai, t'22f = ikN[wi, W2, ^21, ^22 
with N in the form (15. ip . where now 

^2\2 I „2\2 I Q\2 



(6.4) 



T] = cicixi + sixi + 3xi) - 2^2, 

u = C{clXl + slXl)+CXl{l-a2r, 



721 



cx' 
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W2 = a2/{CXl). 



(6.5) 

A search for partial- mode solutions in the form [wi, W2, ^21 5 ^22]"^ = C^"'^*^^, 
where ^ is a 4-component constant vector, leads to an eigenvalue problem. 
The associated characteristic equation is the propagation condition 

{s'-l) [CXls' - C{clXl + slXl) + p] = 0, (6.6) 

with roots s = — 1, which is discarded because it does not lead to a decaying 
wave, 

' (6.7) 



s, = 1, S2 = ±\ [C{clXl + slXl) - p]/{CX 



21 



The analysis leading to the derivation of the dispersion equation is by and 
large the same as that conducted for the principal wave in Section 5. The 
end result is that the dispersion equation is again (15.131) . where now si and 
S2 are given by the expressions above, and 721 is replaced by CA^, but the 
other quantities remain unchanged. 

As an illustration, we take the same numerical values as in Section 5.4, but 
with C2 = in fl5.16l) . Hence the solid is neo-Hookean with C = 0.175 MPa 
and p = 1000 kg/m^. It is under the pre-stress (T2 = 20kPa, (J3 = 0. The 
fluid has Newtonian viscosity fi* = 3.5 x 10~^Ns/m^, and mass density p* = 
1050 kg/m^. The frequency is = 10^ Hz. In turn, we take the solid to be 
compressed by 20% in the Xi direction (so that Ai = 0.8 and then A2 — 1.131, 
A3 ~ 1.105), to be unstretched in the Xi direction (so that Ai = 1, and 
A2 — 1.014, A3 ~ 0.986), and to be under an extension of 20% in the xi 
direction (so that Ai = 1.2, and A2 ~ 0.929, A3 ~ 0.897). 

Figure 4 shows the dependence of the interfacial wave speed = ^{u/k) 
on the angle of propagation 6. When the solid is almost unstrained, the 
speed hardly varies with the angle; when it is strained to ±20%, the induced 
anisotropy causes speed changes of more than ±25% in some directions. The 
figure also shows clearly that the wave travels at its fastest in the direction 
of greatest stretch, and at its slowest in the direction of greatest compres- 
sion, indicating that it could be used for the acoustic determination of these 
directions. 
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